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Abstract 

The disorder and layer transitions in the interface between an Ising spin-1/2 film 
denoted (n), and an Ising spin-1 film denoted (m), are studied using Monte Carlo simula- 
tions. The effects of both an external magnetic field, acting only on the spin-1/2 film , and 
a crystal magnetic field acting only on the spin-1 film, are studied for a fixed temperature 
and selected values of the coupling constant J p between the two films. It is found that for 
large values of the constant J p , the layers of the film (n), as well as those of the film (m), 
undergo a first order layering transition. On the other hand, the only disordered layer 
of the film (n) is that one belonging to the interface films (n)/(m), for any values of the 
crystal field A. We show the existence of a critical value of the crystal field A c , above 
which this particular layer of the film (n) is disordered. We found that A c depends on 
the values of the constant coupling J p between the two films. 
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1 Introduction 

Layered structures consisting of alternating magnetic materials, have been subject of many 
experimental study. These experimental studies have shown that the magnetisation enhance- 
ment exists in multilayered films consisting of magnetic layers. The most commonly studied 
magnetic multilayers are those of ferromagnetic transition metal such as Co or Ni. Several 
experiments have shown that the magnetisation enhancement exists in multilayered films con- 
sisting of magnetic layers. It was found that ferromagnetic coupling can exist between magnetic 
layers. From the theoretical point of view, great interest has been paid to spin wave excitations 
as well as critical phenomena [1-3]. The study of thin films is partly motivated by the devel- 
opment of new growth and characterisation techniques, but perhaps more so by the discovery 
of many exciting new properties, some quite unanticipated. These include, more recently, the 
discovery of enormous values of magnetoresistance in magnetic multilayers far exceeding those 
found in single layer films and the discovery of oscillatory interlayer coupling in transition metal 
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multilayers. 

These experimental studies have motivated much theoretical work. Indeed, several methods 
have been used to study the layering transitions in Ising magnetic films. Benyoussef and Ez- 
Zahraouy have studied these layering transitions using a real space renormalization group [4], 
and transfer matrix methods [5]. Using the mean field theory, Hong [6], have found that depend- 
ing on the values of the exchange integrals near the surface region, the film critical temperature 
may be lower, higher than, or equal to that of the bulk. Using the perturbative theory, Harris 
[7] have showed the existence of layering transitions at T = in the presence of a transverse 
magnetic field. The effect of finite size on such transitions has been studied, in a thin film 
confined between parallel plates or walls by Bruno etal. [8] taking into account the capillary 
condensation effect. By applying Monte Carlo simulations on thin Ising films with competing 
walls, Binder et al. [9], found that occurring phase transitions belong to the universality class 
of the two-dimensional Ising model and found that the transition is shifted to a temperature 
just below the wetting transition of a semi-infinite fluid [10,11]. Hanke etal. [12] showed that 
symmetry breaking fields give rise to nontrivial and long-ranged order parameter profiles for 
critical systems such as fluids, alloys, or magnets confined to wedges. 

We showed in one of our earlier works [13] the existence of layering transitions under the effect 
of a variable surface coupling. Moreover, for an Ising film with a wedge, we found in Ref. [14] 
the intra-layering transitions under the geometry effect consisting on the existence a wedge. 
When the film is subject to a random transverse magnetic field [15], we found the layer-by- 
layer transitions when increasing the concentration p above a critical value p c {k) for each layer 
k. The aim of this work is to study the interface disorder coupled to layering transitions in a 
spin— 1/2 film in interaction with a spin— 1 film, using Monte Carlo simulations. The paper is 
organised as follows. In section 2, we give the model and the method used. In section 3 we 
present results and discussions. 

2 Model and Monte Carlo simulations 

We are studying two coupled magnetic films: a spin-1/2 film denoted (n) formed with n square 
ferromagnetic layers; and a spin-1 film denoted (m) formed with m square ferromagnetic layers, 
Fig. 1. Each layer is a square of dimension N x x N y = 64 x 64 spins. N x and N y being the 
number of spins in the x and y directions, respectively. The juxtaposed plans of the two films 
are coupled ferromagnetically via a constant J p , so that the Hamiltonian governing the system 
is given by: 

H = -J E SiSj-H^Si-J E Wi-Jp E S i a 3 + A E (!) 

<i,j>(n) ie(n) <i,3>(m) ie(n),je{m) je(m) 

where, Si(l = = —1, +1 are the spin variables in the film (n) and a^(k = i,j) = 0, ±1 are 
the spin variables in the film (m). The spins of the film (n) are subject to an external magnetic 
field H, whereas A is the crystal field acting only on spins of the film (m). The sum J2<ij> (n) 
(resp. J2<i,j> (m) ) is performed on nearest neighbour spins of the film (n) (resp. of the film (m)). 
The interaction J between the spins of the film (n) as well as between the spins of the film (m), 
is assumed to be constant. Values of the interface coupling constant J p will be discussed in 
all the following. A preliminary study showed that, when performing Monte Carlo simulations 
under the Metropolis algorithm, the relevant calculated quantities did not change appreciably 
for the film thickness: n,m = 3,4,5 and 8 layers; and when varying the number of spins of 
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each layer from N x — N y — 32 to 128. Where N x and N y are the number of spins, of each 
layer, in the x and y— directions, respectively. Taking into account the above considerations, 
in all the following we will give numerical results for a film thickness n = to = 5 layers, and 
N x = N y = 64 spins for each layer. 



3 Results and discussion 

The geometry of the system we are studying is summarised in Fig. 1, where the two films 
are coupled via the constant J p . Each layer of the film (n) is subject to the external magnetic 
field H; whereas the layers of the film (to) are under the crystal field A. Although the model 
and established equations are valid for arbitrary values number of layers n and to, numerical 
results will be given in this work for n = m = 5 layers. Also the temperature fluctuations will 
not be taken into account and will be fixed at the value T = 3.0 in all the following. 
In order to outline the effect of the external magnetic field, we plot in Figs. 2a and 26 the 
layer magnetisation behaviour in each film: n and to for a fixed crystal field A = 2.0. The 
former figure, plotted for a small coupling constant value J p = 0.5, shows that only the layers of 
film (n) (spin— 1/2) transit when increasing the external magnetic field, including the interface 
layer n(l). Hence for small values of the parameter J p the layering transitions are absent in 
the film (to) (spin— 1) for all values of the external magnetic field. But for higher values of the 
coupling constant J p the layering transitions begin to occur in the film (to), as it is shown in 
Fig. 2b for J p = 9.0. When increasing the amplitude of J p more and more all the layers of the 
film (to) transit, following the layering transitions present in the film (n). This means that for 
sufficiently large values of the coupling constant J P) the effect of the magnetic field H acting 
on the film (n) is propagated to the film (to) layers. 

On the other hand, one can ask if the crystal field A, applied on the film (to) and disordering 
its layers, can affect the layers of the film n. To answer this question, we plot in Figs. 3a 
and 36 the layer magnetisation behaviour in each film as a function of the crystal field A for a 
fixed value of the external magnetic field H = 0.20 and two values of the coupling constant J p . 
Indeed, the only disordering layer of the film (n) is that one belonging to the interface: n(5). 
The disorder of this specific layer is affected neither by small values of the coupling constant: 
J p = 0.5 in Fig. 3a; nor by higher values: J p = 9.0 in Fig. 36. Hence, the disorder of the film 
(to) does not affect the layers of the film (n) except the special layer directly connected to the 
interface between the two films. 

On the other hand, the effect of the interface coupling constant is summarised in Figs. 4a and 
46 corresponding to a low crystal field value A = 0.5 and a higher value A = 9.0, respectively. 
It is found that starting from a totally disordered film (to) for low values of A, when increasing 
the coupling constant J p values, all the layers of the film (to) undergo non null magnetisations 
(keeping the sign of the film (n) layers), see Fig. 4a. While for higher values of the crystal 
field A, see Fig. 46 , the increasing coupling J p values does not affect, neither disorder of the 
film (to) layers, nor the layer magnetisations of the film (n). Indeed, the special layer n(5) 
disordered for the high crystal field value A = 9.0 at low values of the coupling J p , is not 
affected by increasing the coupling J p between the two films. 

In the following we will focus our interest on the particular layers of the interface of the films: 
n(5) belonging to the spin-1/2 film (n), and to(1) belonging to the spin-1 film (to). Indeed, 
concerning the transition of the layers n(5) for a fixed crystal field value A = 1.0 we show in 



3 



Fig. 5a the effect of increasing coupling J p on this layer. It is shown that the transition of the 
layer n(5) inside the film (n) is not affected by increasing the coupling J p from low values 0.1, 
1.0 to a higher value 9.0. Whereas, under the same conditions, the transition of the layer m(l) 
is strongly affected by increasing the coupling J p , See Fig. 56. This layer does not transit for 
J p = 0.1 and any value of the external magnetic field H . 

On the other hand, the layer n(5) disorders for A > 9.0 , see Fig. 6a, for any coupling J p 
value; whereas the layer m(l) disorders for J p = 0.1 at A pa 2.0, for J p = 1.0 at A w 6.0 and 
for J p = 9.0 at A > 10.0, Fig. 66. This means that the crystal field needed to disorder these 
layers depends on the coupling J p of the interface between the two films. A scenario of this 
dependence is plotted in Fig. 7 for the layer n(5) and a fixed external magnetic field H = 0.20. 
It is seen, from this figure, that to disorder the layer n(5) for low values of the coupling J p , 
small values of the crystal field are needed. When increasing the coupling J p , the crystal field 
A c undergoes a maximum value at ~ 9.0, and decreases when increasing the coupling J p more 
and more. 

4 Conclusion 

The disorder and transitions of The interface between an Ising spin-1/2 film denoted (n), and an 
Ising spin-1 film denoted (to), has been studied using Monte Carlo simulations. The effects of 
both an external magnetic field, responsible on the interface transition, and a crystal magnetic 
field needed to disorder the interface are studied for a fixed temperature. It is found that for 
low values of the coupling J p , the layers of the film (n) undergo a first order layering transition. 
These transitions are also found in the film (to), for higher values of the coupling J p between 
the two films. On the other hand, the layers of the film (to) disorder when increasing the crystal 
magnetic field for any coupling J p value, but the only disordering layer of the film (n) is that 
one belonging to the interface between the films (n)/(m). 

On the other hand, we show the existence of a critical value of the crystal field A c , above which 
this particular layer of the film (n) is disordered at given values of the interface coupling J p . 
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Figure Captions 
Figure 1.: 

Geometry of the system formed with two films coupled via the constant J p . Each layer of the 
film (n) is subject the external magnetic field H; whereas the crystal field A is acting only on the 
layers of the film (m). The films (n) and (m) are formed with the same number of layers, N = 5. 

Figure 2: 

The external magnetic field effect on the layer magnetisation in each film, for a fixed crystal 
field A = 2.0. (a) For a small coupling constant value J p = 0.5, only the layers of film (n) 
transit when increasing the external magnetic field, including the interface layer n(5). The 
layering transitions are absent in the film (m) for this small value of the coupling constant, (b) 
For a higher value of the coupling constant J p = 9.0 the layering transitions begin to occur in 
the film (m) following the layering transitions present in the film (n). 

Figure 3: 

The layer magnetisation behaviour in each film as a function of the crystal field A for H = 0.20. 
(a) For a small value of the coupling constant: J p = 0.5, when increasing A all the layers of 
the film (m) are disordered but the only disordering layer of the film (n) is that one belonging 
to the interface: n(5). (b) For a higher value of the coupling constant: J p = 9.0, the disorder 
of the film (m) layers does not affect the layers of the film (n) except the special layer n(5) 
directly connected to the interface between the two films. 

Figure 4: 

The interface coupling effect on the layer magnetisation behaviour in each film for a constant 
external magnetic field H = 0.20 and two crystal field values: (a) For A = 0.5 and starting from 
a totally disordered film (m) all the layers of the film (m) undergo non null magnetisations when 
the coupling J p between the films is increasing, (b) For a higher crystal field value A = 9.0, the 
increasing coupling J p values does not affect, neither the disorder of the film (m) layers, nor 
the layer magnetisations of the film (n), except the special layer n(5) belonging to the interface. 

Figure 5: 

Magnetisation behaviour of the interface layers: n(5) and m(l) for a fixed crystal field value 
A = 1.0 and different values of the coupling constant J p = 0.1, 1.0 and 9.0, as a function of 
the external magnetic field H. (a) The transition of the layer n(5) inside the film (n) is not 
affected by increasing the coupling constant, (b) The layer m(l) does not transit for low values 
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of the coupling, e.g. J p = 0.1. The transition of this layer is only seen for large values of the 
coupling J p . 

Figure 6: 

Magnetisation behaviour of the interface layers: n(5) and m(l) for a fixed external magnetic 
field value H = 0.20 and different values of the coupling constant J p = 0.1, 1.0 and 9.0, as a 
function of the crystal field A. The crystal field needed to disorder the layer n(5) (a) and the 
layer m(l) (b) depends strongly on the coupling J p between the two films. 

Figure 7: 

The critical crystal field A c needed to disorder the layer n(5) as a function of the coupling J p , 
for a fixed external magnetic field H = 0.20. A c presents a maximum value for a coupling 
J p w 12.0 and decreases when increasing the coupling J p between the two films. 
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